Spectroscopy of a Cooper-Pair box in the Autler-Townes configuration 
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Abstract: A theoretical spectroscopic analysis of a microwave driven superconducting charge 
qubit (Cooper pair box) coupled to an LC oscillator model is performed. By treating the oscillator 
as a probe through the backreaction effect of the qubit on the oscillator circuit, we extract frequency 
splitting features analogous to the Autler-Townes effect from quantum optics, thereby extending the 
analogies between superconducting and quantum optical phenomenology. These features are found 
in a frequency band that avoids the need for high frequency measurement systems and therefore may 
be of use in qubit characterisation and novel coupling schemes. In addition we find this frequency 
band can be adjusted to suit an experimental frequency regime by changing the oscillator frequency. 
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I. INTRODUCTION 

Superconducting charge qubits (Cooper pair boxes) 
are promising candidates for use in quantum informa- 
tion processing and quantum computing systems. Many 
of the necessary results for quantum computing have al- 
ready been demonstrated in these systems including sin- 
gle qubit rotations [10] , single-shot readout and two- 
qubit entanglement jjj. 

In common with other solid-state technologies (e.g. 
Si:P superconducting devices are expected to be 

highly scalable due to the ability to construct large num- 
bers of qubits using existing lithographic manufactur- 
ing 0. However, any lithographic procedure has in- 
trinsic manufacturing tolerances and there can be sig- 
nificant variation in device parameters between neigh- 
bouring qubits. These are of little consequence in classi- 
cal digital logic designs because the 'zero' and 'one' logi- 
cal states are usually well separated in terms of the lev- 
els used to represent them 0. They are well-defined 
discrimination thresholds and they are highly robust to 
small fluctuations. Such classical error protection fails 
as the number of particles (electrons) defining a given 
state is reduced to the few particle level and the fragility 
of quantum coherence introduces new problems. Hence 
the Hamiltonian of each individual qubit will need to be 
characterized to unprecedented levels 0, Q . 

Recently, attention has turned to the coupling of su- 
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perconducting flux and charge qubits to microwave os- 
cillators and resonators [13 111 III 111 13 These will 
have some resonant frequencies caused by the intrinsic 
capacitance and inductances in the connection path, this 
can often be approximated by linear oscillators for small 
fluctuations and deviations. The simplest example of this 
is in the usual spectroscopy of the qubit where the mi- 
crowave field acts as a weak probe |l5), and increasing 
the strength of the field allows Rabi oscillations 0] . To 
investigate the strong-coupling limit of cavity Quantum 
Electrodynamics, Wallraff et al. constructed the analog 
of an optical cavity for a Cooper-Pair Box, and showed 
single photon dynamics |T^ . This work and related the- 
oretical analysis Hal c learly suggests that the full canon 
of quantum optics pjj can be translated to the supercon- 
ducting arena for potential technological advantage. In 
addition, the ability to fabricate desired pseudo-atoms 
allows investigation of phenomena (particularly multi- 
state phenomena) that would not be convenient using 
real atoms. 

We have investigated a simulated superconducting 
charge qubit model, controlled via a classical microwave 
pump and a bias/control voltage. We treat the bias con- 
trol as a fully quantized linear oscillator with a finite re- 
sistance (RLC oscillator) [2(j which is coupled to a ther- 
mal environment |2l| . The quantised oscillator is used to 
allow the field to represent non-adiabatic processes and 
entanglement between the qubit and bias field. 

System measurements can be performed on the oscilla- 
tor circuit or the qubit, but measurement of the oscillator 
is generally preferred as it is more readily accessible via 
weak measurements than the fragile qubit state. One of 



2 



osc 




QUBIT 



Oscillator Qubit 



mw 



= 3.49GHz 



mw 



= 3.49GHz 



Fnsc — 3.06GHz 



Low 
frequency 
oscillator 



Fnsc= 1.36GHz 



High 
frequency 
oscillator 



FIG. 1: A two level qubit is coupled to a many level harmonic oscillator, investigated for two different oscillator energies. 
Firstly, the oscillator resonant frequency is set to 1.36GHz, this more resembles the conventional configuration such that the 
fundamental component of the oscillator does not drive the qubit. However, we also investigate the use of a high frequency 
oscillator of 3.06GHz which can excite this qubit. In addition, qubit is constantly driven by a microwave field at 3.49GHz to 
generate Rabi oscillations and in this paper we examine the relation between these three fields. 



the ways to reduce decoherence is to reduce the num- 
ber of degrees of freedom coupling directly to the qubit, 
thereby limiting the environmental noise. In addition, 
the oscillator is of particular interest as its evolution is 
significantly affected by the qubit dynamics due to back- 
action caused by the two way capacitive qubit-oscillator 
coupling. 

This paper examines the power spectral density (PSD 
or noise spectrum) of the oscillator and qubit as functions 
of the two main control parameters, microwave amplitude 
and bias voltage. The PSD, particularly for low sub GHz 
frequencies, is a readily accessible experimental measure- 
ment |22l l23j| . By exploring non-standard parameters for 
the bias circuitry and microwave control, our investiga- 
tions show a rich spectroscopy of dynamically control- 
lable features. Beyond fundamental investigations, these 
results may be of benefit in the characterisation of qubit 
parameters, and the tunability of the features may al- 
low for new control techniques for single and multiple 
qubit coupling 24]. These features are not restricted to 
superconductive devices, as semiconductor charge qubits 
may also exhibit similar behaviour. The magnitude of 
the features are of a smiliar order to the oscillator peak 
and so should be detectable if the oscillator peak can be 
sucessfully measured. 

We consider two situations (Fig.[TJ in which we change 
the resonant frequency of the biasing circuit, we have cho- 
sen the values to be non-commensurate as to avoid the 



possiblity of higher oscillator harmonics repeating the 
same effects. Initially, we examine a conventional con- 
figuration in which the bias oscillator frequency (/ osc = 
1.36GHz) is much less than both the undriven qubit tran- 
sition frequency (f qu bit = 1.94GHz) and the microwave 
drive frequency (f mw — 3.49GHz). This oscillator fre- 
quency is too slow to create excitations, so only a single 
set of microwave driven excitations (Rabi oscillations) are 
observed. The second scenario, involves using a high fre- 
quency oscillator (f osc — 3.06GHz) (2j| near the qubit 
and microwave frequencies. This smaller detuning cre- 
ates a range of interesting features in the oscillator spec- 
trum. We also observe a frequency splitting when the 
qubit is correctly biased in the presence of a driving signal 
and strong coupling. This resembles the Autler-Townes 
effect observed in quantum optics |26j. In this case, the 
oscillator frequency is sufficient to create its own exci- 
tations, which intersects and mixes with the microwave 
driven Rabi oscillations. This mixing allows the Rabi fre- 
quency to be modulated/upconverted into an experimen- 
tally favourable frequency band, in addition we find that 
the position of this band is can be adjusted by changing 
the oscillator frequency. 

The Autler-Townes (A-T) experiment is a well-known 
experimental scheme from quantum optics for probing 
a strongly driven transition. In the canonical example 
|20j . an atomic transition is driven by a strong, quasi- 
resonant field. This induces Rabi oscillations, and due 
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to decoherence processes, and if the driving is resonant, 
the population in the two states will tend to equalise. In 
this case the most usual way to describe the system is 
in the dressed-state basis, where the dressed states are 
the eigenstates of the driven system, and correspond to 
orthogonal coherent superpositions of ground and excited 
state. The transient evolution of a three-state system in 
an A-T-style configuration is studied experimentally and 
theoretically in [271 ]. 

Because the population in ground and excited state 
tends to equalise in the long time limit, the fidelity 
of a readout signal can be compromised. Against this 
backdrop, the A-T configuration provides an additional 
readout handle. By sweeping a weak probe field, quasi- 
resonant with a transition from the one state to a third, 
otherwise uncoupled state, and ensuring that the sweep 
range covers all of the substructure of the dressed states 
induced by the strong coupling field, the population in 
the dressed states, rather than the bare atomic states is 
read out. As this is a more natural basis to be monitor- 
ing in a strongly driven system anyway, the experiment 
allows a useful measurement of otherwise hidden dynam- 
ics to be performed. Good examples of systems well 
resolved by this process are bichromatic |23 and poly- 
chromatically driven transitions [2^, |3(| , and V-systems 
[3l| . In general because the third state is assumed un- 
affected by the strong driving, then the population dif- 
ference between the dressed states and the third state 
can, in principle, be very large, and is unaffected by the 
strong coupling. 

More generally, however, the A-T experiment provides 
for more flexibility, and introduces an alternative mea- 
surement handle. In this context it has potential utility 
in solid-state systems, where the Rabi oscillations are 
replaced with coherent-tunnelling oscillations, and the 
tuning of transitions is performed electrically via sur- 
face gates. In this paradigm, schemes related to the A-T 
experiment have been proposed as alternative readout 
mechanisms for a charge qubit operating in the superpo- 
sition basis |32j | and for singlet-triplet discrimination in 
a two-spin system 33]. As will be shown below, the in- 
creased flexibility of the A-T configuration in our system 
emphasises side-bands corresponding to subharmonic res- 
onances and features at multiples of the Rabi fre- 
quency frequency [33 , which would otherwise be hidden 
in more conventional, two-state systems, and afford new 
diagnostic methods, and possibilities for control. 



II. SYSTEM MODEL 

The system model (shown in Fig. [2J is a small su- 
perconducting island, capacitively coupled to a bulk su- 
perconductor via a Josephson junction, allowing Cooper 
pairs to coherently tunnel on and off the island. The 
Josephson junction is sufficiently weak and the total ca- 
pacitance of the island sufficiently small that the energy 
associated with charging the island is large compared to 



the tunnelling energy of the junction (hv\ this is nor- 
mally called the Coulomb blockade regime 36] . The nec- 
essary biasing voltage is applied through the bulk super- 
conductor, across the Josephson junction. In addition, 
the qubit is also coupled capacitively to a grounding elec- 
trode, a reference point for the junction biasing voltage. 

The bias circuit has natural intrinsic dynamics with 
a characteristic frequency response, we approximate this 
by a resonant RLC oscillator in series with the island, 
which takes into account the intrinsic inductance and ca- 
pacitance present in the bias circuit, as well as a par- 
asitic capacitance across the bias electrodes. This os- 
cillator is modelled quantum mechanically as a standard 
TV-state Simple Harmonic Oscillator (SHO) allowing non- 
adiabatic processes and not limiting the analysis to low 
frequency fields |37ll38tl39tl4(| . The bias oscillator is also 
connected to an external amplifier through which the os- 
cillator voltage can eventually be measured. We have 
focused on the effects of components immediately cou- 
pling to the qubit, although it should be acknowledged 
that the specific entire chain 41] which will add more re- 
sistive, capacitive and inductive loads. The qubit can be 
indirectly measured using its reactive backreaction on the 
oscillator voltage, the backreaction is a result of the two- 
way capacitive coupling between the qubit and oscillator 
circuit. Indeed, the presence of fluctuations in the oscil- 
lator will also have an effect on the qubit, an additional 
source of environmental decoherence, this is minimised 
by maintaining a suitably low temperature (lOmK). 

The Hamiltonian for this coupled system is: 



H = —— + K— + —— + V in Q 
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with Q and q being the effective oscillator and island 
charges respectively, $ is the inductor flux and 8 the su- 
perconducting phase difference, the effective island and 
bias electrode capacitances and coupling parameter k are 
given by: 
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The qubit geometry is transformed to a classical elec- 
trical circuit model (Fig. |2J). The circuit includes the 
intrinsic qubit capacitances Cj and C g , (the Josephson 
junction capacitance and a grounded bulk connection ca- 
pacitance). In addition, a large parasitic capacitance C p 
is formed between the bulk superconductor electrodes 
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FIG. 2: An abstract layout of a simple single junction Cooper pair box, with the biasing voltage connection to the first filter 
stage modelled by a resistor and inductor. The Cooper pair box is modeled by an equivalent circuit network, consisting of three 
capacitors: C.j - the capacitance of the Josephson junction, C g - The capacitance of the island gate, and C p - the parasitic 
capacitance formed between the biasing electrodes. These capacitors, together with the resistor and inductor form an RLC 
oscillator circuit, the inductance value is changed to investigate the two oscillator frequencies, (1.36GHz and 3.06GHz). 



TABLE I: Components (Capacitors same for both fosc) 





Description 


Value 


v/2-k 


Josephson energy 


1.94GHz 


Cj 


Josephson junction capacitance 


0.5fF 




Qubit-Grounded Bulk capacitance 


0.5fF 


Cp 


Electrodes parasitic capacitance 


lOpF 


Lose 


Bias inductance (f OS c = 1.36GHz) 


1.37nH 


Rose 


Bias resistance (fosc = 1.36GHz) 


60mft 


7 


Qubit jump rate {fosc = 1.36GHz) 


6.8 x 10 7 


Tl 


Relaxation time (fosc = 1.36GHz) 


14.7ns 


Lose 


Bias inductance (fosc = 3.06GHz) 


0.27nH 


Rose 


Bias resistance (fosc = 3.06GHz) 


27mft 


7 


Qubit jump rate (fosc = 3.06GHz) 


1.53 x 10 8 


Tl 


Relaxation time (fosc = 3.06GHz) 


6.5ns 


K 


Coupling parameter 


5 x 10" 5 


Qd 


Oscillator quality factor 


200 



and is connected to the inductive and resistive elements 
of the biasing circuit. 

The Hamiltonian is derived from consideration of cir- 
cuit voltage and current laws using a detailed loop analy- 
sis method proposed by Burkard [42] , with the Josephson 
junction current given by the well known phase-current 
relation y| where 9 is the superconducting phase differ- 
ence between the adjoining segments. 

/, = / c sin( 27 r£) (3) 

The resulting Hamiltonian \Q is rewritten to separate 
the qubit (q) and oscillator (Q) degrees of freedom, and 
a coupling constant k is introduced. The k term is typ- 
ically much less than one, and in the limit k — > the 
cross-coupling qQ term is zero and the systems evolve 
independently. 



In addition, it should be noted that the capacitances 
Cj, C g , C p have been redefined as effective qubit and bias 
oscillator capacitances, C q , C q Q,CQ, ||2J). The term Vi n Q 
is the contribution of the bias/control voltage and the 
inductive term arises from the magnetic flux associated 
with the oscillator geometry. 

Starting from Eq. Q it is possible to group the Hamil- 
tonian terms into three parts: the oscillator and qubit 
Hamiltonians, plus a coupling term. 

H H osc C*D Iqubit lose ® Hqubit ~t~ H coupling (4) 

where I is the corresponding identity matrix for each 
component. The Hamiltonian is written in the matrix 
representation, using a basis formed from the tensor 
product of the uncoupled oscillator energy eigenstates 
and the charge states of the isolated qubit (which 
forms the compuational basis for the qubit). The three 
terms in the Hamiltonian are presented below in Eqs. 
@ El nT)l . In the these equations, the subscripts are 
used to label the basis states of the uncoupled systems: 
the TV-state oscillator is indicated by a subscript T' and 
the qubit by a subscript '2'. So that m\ or n\ represent 
the oscillator states (mi,ni € {0, 1,2, ...,7V}) and 7712 or 
712 indicate the qubit states (111,2, ^2 € {0, 1}). 

1. Oscillator Hamiltonian, H osc 

The oscillator is quantised as a simple harmonic 
oscillator JSJ ^2|> with energy oj. In the electri- 
cal circuit model, the charge Q and flux $ form 
the conjugate variables, analogous to momentum 
(p) and position (x), (|5a.p. The usual mechanical 
model of the harmonic oscillator can be equated 
with the electrical Hamiltonian, yielding equivalent 
values for an effective 'mass' m = Cq and 'spring 
constant' k = l/L osc . We also assume if the qubit is 
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maintained at a suitably low temperature (lOmK) 
the contributions of the higher energy levels are 
negligible, hence we constrain the system to the 
lowest N = 10 oscillator states for computational 
efficiency, whilst checking that probabilities are not 
lost to the higher states. 



Hose 


$ 2 




2Losc 






= Hui n -\ 






tensoring with the qubit basis yields Eq. ©, the 
first part of the combined system matrix Q 



(m 1 ,m 2 \ H osc \n 1 ,n 2 ) 
fku \n\ + — 



^mi ,ni ^wt2 ,ri2 (^) 



where, in this paper, u/2'k is 1.36GHz for the low 
frequency oscillator with 3.06GHz being the high 
oscillator frequency. 

2. Qubit Hamiltonian, H qu bit 

The qubit basis is defined as the presence or ab- 
sence of a single Cooper pair on the island, Oe and 
2e. Although there may be higher states present, 
we assume the probability of a third state is negli- 
gible 



H, 



qubit 



hqubit 



2C a 



hvcos 2ir — 



|L|i)(i|-^(| ) (i| 



ID 



(7a) 
(7b) 



where the cosine term can be written as the sum of 
two exponentials divided by two, these exponentials 
are raising and lowering operators on charge states, 
which are then expressed in the matrix formulation 
as off diagonal terms. Tensoring this small matrix 
with the oscillator basis yields Eq. ©, the second 
part of the combined system matrix Q . 



(mi, 771 2 | H qubit |t7i,77 2 ) 



(8) 



The biasing field caused by voltage Vi n can be rep- 
resented as an effective bias charge Quas ■ The mi- 
crowave field is applied as an effective qubit charge, 
the combination is q ext 0]. 

sin (oj mw t) (9a) 
q(t) = gext|0)<0| + (ft nrt + 2e)|l)<l| (9b) 



3. Coupling term, H coup i ing 

The qubit-oscillator coupling term is the product 
of the qubit and oscillator charge operators. To 
quantise this we express the oscillator charge Q in 
terms of raising and lowering operators 20]. The 
terms 



fn{ and \Jn\ + 1 are normalisations for the 
changes in oscillator state. 



H, 



coupling 



Q 



qQ 



hmuj 



(10a) 
(10b) 



the tensor product of the q and Q matrices yields 
Eq. (|lll) . the final part of the combined system 
matrix Q). 



(m 1 ,m 2 \ H qubit \n 1 ,n 2 ) 



C 



\j2fvmu) 8 m2 iS„ 



III. 



(y/n\ 6m un -i - Vni + 1 5 mi , n +i) (H) 



UNRAVELLINGS AND WEAK 
MEASUREMENTS 



The evolution of the quantum system (qubit and os- 
cillator) in the presence of environmental fluctuations is 
described by the usual master equation. 



dp 



-i[H,p) 

( L ™P L tn - \ L \n L mP- \pL nl L m \{\2) 



where the ensemble average of the individual state evolu- 
tions yield the density operator for the combined system 
(The qubit and oscillator circuit) p = E (\ip) (ip\). Eq. 
(|12[) is the general equation for the average evolution of 
an ensemble ^i), for which the Caldeira-Leggett model 
[i^j l is a specific example. The unravellings represent the 
evolution of specific systems (not ensemble evolutions) 
under the actions of projective measurements on their 
environment. Here we have ignored the effects of expicit 
decoherence due to dephasing in individual qubits. 

Unravelling the Master equation allows us to represent 
only the stochastic effects within the qubit and oscilla- 
tor as a whole. In this section we describe two unrav- 
ellings which correspond to spontaneous emissions and 
the thermal diffusion mechanism of the oscillator. This 
environment also includes the effects of an external mea- 
suring system on the oscillator as this will add to the 
decoherence. We have chosen a quantum state diffusion 
(QSD) unravelling (|13|) to model oscillator dissipation, 
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although equally valid is the quantum jump (QJ) unrav- 
elling representing a different measurement process on 
the environment, indeed both unravellings should agree 
on average 0. 



A. Oscillator quantum state diffusion 

Quantum state diffusion is applied to simulate the ef- 
fects of a finite temperature environment, and dissipation 
|2l| . The diffusion model consists of random thermal and 
quantum fluctuations combined with a progressive drift 
of the quantum state. The system stochastic evolution 
is now governed by Eq. Ijl^JI . The first term is the usual 
Schrodinger evolution, the second term is deterministic 
and represents the drift, and the final term models the 
random fluctuations of the state. The intensity of these 
effects is governed by the Lindblad operators defined in 

Eq. inni 



\chl>) 



-H lib) dt 



J2 (<4^ L ™ - \tf m L m ~ (Ll)^ (L m )\ |^) dt 

■J2(L m -{L m ))\ili)dZ m (13) 



where d^ m are independent complex differential random 
variables (Weiner increments), which satisfy the following 
statistical averages: 



E (dU) 
E(d£ n d£ m ) 

E(dCdU) 

For the arbitrary index m = 
ment operators L m — L\,L,2, 








(14a) 
(14b) 
(14c) 



1, 2 we define two environ- 



Li 
L 2 



(n + 1) — a 



n Q- D a 



exp(^) 



1 



(15a) 
(15b) 
(15c) 



where Qd is the quality factor of the oscillator and is 
proportional to the dissipative elements, n is the aver- 
age oscillator occupancy and lu is the angular frequency 
of the oscillator. The oscillator occupancy is dependent 
on the thermal temperature of the environment which is 
assumed constant at T=10mK. 



B. Qubit quantum jumps 

To model the spontaneous photon emissions from the 
qubit as it undergoes a modified Schrodinger evolution 



we adapt a quantum jump trajectory unravelling found 
in quantum optics |47|]. We allow the qubit to emit a 
photon into a dissipative reservoir (experimental cavity), 
in which it is immediately absorbed (An ideal measure- 
ment, to simulate a quantum jump). If we apply con- 
tinuous quantum measurements to the qubit we require 
two measurement operators ()16|l. the subscript indicates 
the emission and absorption of a photon (i.e. a jump has 
occurred). 



Q (dt) = 1 - ^Hdt - ]-<r^adt 
fli (dt) = \J jdt a 



(16a) 
(16b) 



where a = |0) (1|, while 7 sets the jump rate so that 
7 (<rV) is the actual photon emission rate, (7 = 0.05/ osc 
which implies 0.05 jumps per bias oscillator period, note 
that 7 differs for the two bias oscillator frequencies, see 
Table QJ. Eq. (|16a|l represents the continuous evolution 
of the qubit state over the interval dt, whereas Eq. (|16b|) 
describes the collapse of the state when a quantum jump 
has occurred. Over each simulation interval, dt, the state 
is updated via one of two Eqs. H17fl dependent on the 
random result of an instantaneous weak measurement, 
in normal evolution Eq. i|17al) is used, however if a jump 
occurs a '1' will be measured (photon emitted) and Eq. 
(|17b(l should then be used and the state renormalised. 



IV. 



|Vo (* + <&)) = fl (dt)\i/j(t)) 

IVi (* + *)> = ni(dfc)|V(t)> 



ENERGY LEVEL STRUCTURE 



(17a) 
(17b) 



It is useful to be able to quantify several important 
qubit parameters, as the qubit behaviour is particularly 
sensitive to parameter variation. These arise from the 
manufacturing techniques used, which have often been 
driven by industrial goals of delivering cheap yet inac- 
curate components. It is important to understand that 
small variations in the capacitive coupling can have a 
significant effect on the system behaviour. These errors 
could occur through a variety of real situations. For ex- 
ample, overlaps between layers in the epitaxially grown 
structure, intrinsic parasitic capacitances such as fringe 
effects, or defects in the Josephson junctions |4g|. Un- 
fortunately whilst these variations can be kept to a min- 
imum by improved design and lithography, we cannot 
expect perfect structures, and once integrated and fabri- 
cated we cannot measure every capacitance directly. In- 
stead, we are likely to require methods to remotely char- 
acterise the qubit in circuit. 

The qubit behaviour is partially characterised by the 
energy level diagram (Fig. [Hi EH, which shows the 
time averaged (Floquet) energy levels for the microwave 
driven qubit [5(J, |5ll . The Floquet energies are a useful 
illustration of where transitions can occur for a particu- 
lar microwave drive frequency, and show how the qubit 
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0.45 0.475 0.5 0.525 0.55 

Bias, (2e) 

FIG. 3: (Color online) Time averaged Floquet qubit energy levels, with one and two photon transitions (3.49GHz and 6.98GHz), 
following results focus on the region surrounding the single photon transition (3.49GHz) at Bias = 0.5187. The solid upper and 
lower lines are the time averaged energies of the two states during microwave drive, the dashed upper and lower lines are the 
undisturbed energies. The energy diagram clearly indicates the bias values near which the microwave transitions exist. Time 
averaged Floquet energies are necessary to conveniently show the transitions as the microwave drive is a time varying field. 



energy eigenvalues change as a function of the microwave 
and bias control fields. However, this is not necessar- 
ily known and cannot be assumed to be identical for all 
qubits. 

Indeed, techniques have been proposed to deduce the 
energy level structure from measurements of the in- 
creased peak noise power in the biasing circuit, caused by 
Rabi oscillations when the qubit is correctly biased (on 
resonance) with the injected microwave field [371 l38l l39| . 

The system described in this paper is driven with a 
3.49GHz microwave, which gives rise to a transition be- 
tween the energy levels, the width of which becomes 
wider as the microwave amplitude (A mw ) is increased. 
The single photon (3.49GHz) transition occurs symmet- 
rically near two bias points (0.4812 and 0.5187). We have 
chosen to operate about the upper bias point (0.5187), 
although identical features would be present near the 
lower point. Two photon transitions occur at 6.98GHz, 
although the bias required to utilise these transitions are 
far above our parameter range, so their effects will be 
ignored. 

The graph also clearly shows the expected avoided 
crossing caused by the Josephson junction, with an en- 
ergy separation, Ej (1.94GHz), transitions cannot oc- 
cur for microwaves driving below this frequency. This is 
important, as only the 3.06GHz high frequency biasing 



oscillator should be sufficient to create the excitations 
that interfere with the microwave and cause the many 
strong and interesting features, unlike the 1.36GHz low 
frequency oscillator circuit. 

V. FREQUENCY SPECTRA 

In this section we consider the frequency spectrum of 
the system, or the power spectral density (PSD) which 
is an experimentally measurable quantity. The power 
spectral density reveals the strength of the frequencies 
present in a time varying signal. In this case the signals 
are voltages, which are calculated from the expectation 
values of qubit and oscillator charge. The expectation 
values are obtained using the reduced density operator 
for the qubit and oscillator basis, applying the relevant 
charge operator and then tracing the result. 
E.g. (q) = Tr {qp q } where p q is the reduced density 
matrix for the qubit state. 

P(Q) = 201og 10 \FFT(V Q )\ (18a) 
P(q) = 201og 10 \FFT(V q )\ (18b) 

where Vq = (Q) /C Q and V q = (q) /C q . Vq should 
be accessible though measurements on the biasing cir- 
cuit, and V q (if required) can be measured using a Single 
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Oscillator PSD as a function of bias A. ... = 0.0025 Qubit PSD as a function of bias A,.„, = 0.0025 




Frequency (GHz) Frequency (GHz) 

FIG. 4: (Color online) Oscillator and Qubit power spectra slices for Bias = 0.5187, using the low frequency oscillator circuit 
fosc = 1.36GHz. The solid lines overlay the energy level separations found in Fig- 01 = 5 x 1CF 5 ). As one would expect, the 
bias oscillator peak at 1.36GHz is clearly observed in the oscillator PSD, but only weakly in the qubit PSD. Likewise the qubit 
Rabi frequency is found to be stronger in the qubit PSD. However it is important to note that the qubit dynamics such as the 
Rabi oscillations are indeed coupled to the bias oscillator circuit and so can be extracted. In addition, it is recommended to 
compare the layout of the most prominent features with Fig. |S| 



Electron Transistor (SET). FFT denotes the Fast Fourier 
Transform which is a computationally efficient implemen- 
tation of the standard Fourier transform, to obtain the 
frequency spectra of the stochastic time evolutions Vq 
and V q . 

In the example figure (Fig. 0}, the control bias is var- 
ied from left to right for a low frequency oscillator cir- 
cuit (1.36GHz). For each bias point the simulation is 
reinitialised, the stochastic time evolution of the sys- 
tem density matrix is simulated over 1500 oscillator cy- 
cles. Then the oscillator and qubit charge expectation 
values are extracted to obtain the power spectrum for 
each component, with a frequency resolution of 4.01MHz. 
The power spectra for each time series are collated as 
an image such that the power axis is now represented 
as a colour, and the individual power spectra are verti- 



cal 'slices' through the image. The dominant frequency 
peaks become line traces, therefore illustrating the vari- 
ous avoided crossings, mergeings and intersections. The 
example figure shows the PSD 'slice' at Bias = 0.5187, 
the broadband noise is readily apparent and is due to 
the discontinuous quantum jumps in the qubit. The bias 
oscillator peak (1.36GHz) is most prominent in the oscil- 
lator PSD, as would be expected, but it is also present 
in the qubit PSD. It should also be noted that most fea- 
tures are present in both the qubit and oscillator, includ- 
ing the noise which is generated by the quantum jumps 
and the quantum state diffusion processes. Interestingly, 
the qubit PSD is significantly stronger than the oscilla- 
tor PSD, however, a larger voltage is generated by the 
smaller charge due to the extremely small island capaci- 
tance, (V q = q/Cq). 
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Oscillator PSD as a function of bias A. ... = 0.0025 




.51 0.515 0.52 0.525 0.53 0.535 0.54 




Bias, (2e) 

FIG. 5: (Color online) (A) Oscillator power spectra when the coupled qubit is driven at f mw — 5.00GHz. An increase in 
bias noise power (/ osc = 1.36GHz) can be observed when Rabi oscillations occur, the more frequent quantum jump noise 
couples back to the oscillator. (B) Bias noise power peak position changes as a function of fmw, the microwave drive frequency. 
Therefore, it is possible to probe the qubit energy level structure by using the power increase in the oscillator which is already 
in place, eliminating the need for additional measurement devices. However, it should be noted that the surrounding oscillator 
harmonics may mask the microwave driven peak. (« = 1 x 10~ 3 ). 



In addition, the energy level separation (Fig. |3J) has 
been overlaid as a solid curved line to illustrate the effect 
of the microwave drive. 

In a previous paper [37j, a method was proposed by 
which the energy level structure of a charge qubit can 
be obtained from measurements of the peak noise in the 
bias/control oscillator, without the need of extra readout 
devices. This was based on a technique originally pro- 
posed for superconducting flux qubits [37j but there are 
many similarities between the two technologies. The os- 
cillator noise peak is the result of broadband noise caused 
by quantum jumps in the qubit being coupled back to the 
oscillator circuit. This increase in the jump rate becomes 
a maximum when the Rabi oscillations are at peak ampli- 
tude, this should only occur when the qubit is correctly 
biased and the microwave drive is driving at the tran- 
sition frequency. Therefore by monitoring this peak as 
a function of bias, we can associate a bias position with 
a microwave frequency equal to that of the energy gap, 
hence constructing the energy diagram (Fig. |3J). 

The peak is approximately -1-lOdB to +15dB above the 



background noise (Fig. JSJ, which is set by the thermal 
fluctuations. Although it would be difficult to measure 
directly in an experiment, the electronic noise temper- 
ature at the first (low temperature) amplification stage 
should not be significantly above this thermal level or 
it would adversely affect the ability to measure quantum 
mechanical behaviour in the qubit. Care should be taken 
to avoid detecting a fixed peak, these extra peaks are due 
to oscillator harmonics driving the qubit, however the far 
right peak is due to resonance between the Rabi oscilla- 
tions and the bias oscillator, (Fig. 03). 

The next two sections illustrate how the two control 
fields affect the features, by first using the above men- 
tioned method for determining the energy level struc- 
ture, it may be possible to determine the true bias and 
microwave fields applied to the qubit: 

• Section A investigates the effect of sweeping the 
bias field about the energy level transition, for a 
constant microwave frequency of 3.49GHz. 

• Section B investigates the effect of sweeping the 
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microwave frequency, at a constant bias of 0.5187. 

In both sections, the example figures provided use the 
high frequency (3.06GHz) bias oscillator. This higher 
frequency bias field requires a quantised oscillator model 
because the classical model is insufficient |40|. To illus- 
trate the additional effects created when two frequencies 
are driving in close proximity, it is necessary to compare 
the following plots with Fig. 0] In addition the coupling, 
k is kept constant at 5 x 10 -5 . 

A. Bias field sweep 

Fig. shows plots of the PSD for two different mi- 
crowave amplitudes A mw — 0.0025 (A) and A mw = 
0.0050 (B) driving at 3.49GHz with a high frequency bias 
oscillator of 3.06GHz. 

Although not shown, when the microwave field is re- 
moved there are still some effects due to the 3.06GHz 
high frequency oscillator causing excitations and a split- 
ting in the qubit. This is possible as the oscillator ex- 
ceeds the Josephson junction frequency (1.94GHz) and 
so drives the qubit as per an injected microwave field. 
Alternatively, if a low frequency oscillator (1.36GHz) is 
used, there are no excitations and the qubit transition 
frequency remains undisturbed, unless the qubit is bi- 
ased on a harmonic of the oscillator, which would there- 
fore exceed the Josephson frequency and start driving the 
qubit. 

When the microwave is applied we notice interesting 
effects, arising from the additional interference of the 
'side bands' created by the 3.06GHz oscillator drive sig- 
nal, which is now in close proximity to the 3.49GHz 
microwave field. This relatively small detuning causes 
a multitude of splittings, similar to that of the Autler- 
Townes effect [2£|. However, this is a benefit as it creates 
many more features in the frequency spectrum, which in 
turn increases the frequency bands in which useful in- 
formation can be obtained. Although the features ob- 
tained are clearly related, the presence of multiple fre- 
quencies (mixtures, splittings and harmonics) in different 
regions of frequency provides different experimental 'win- 
dows' within which to explore the behaviour of Josephson 
charge devices. 

The mixing of the Rabi oscillations with the microwave 
field causes 'side bands' surrounding the drive signal. It 
is important to note that the high frequency behaviour 
surrounding the drive signals is replicated at f m w — fosc = 
430MHz, and this area of interest has been enlarged, in- 
deed, this should be an easily accessible frequency range. 
This frequency regime can be adjusted on to an experi- 
mental bandwidth by changing the oscillator frequency, 
and hence the aforementioned separation, f mw — f osc . 

In addition, within the PSD there exists a point of 
interest (<0) where the microwave and oscillator excita- 
tions intersect. As the hyperbolic shape of the oscillator 
excitations remains quite constant for low amplitude mi- 
crowaves, the intersection point will mainly be governed 



by the amplitude and frequency of the microwave. It can 
be seen in plot A, when the microwave is driving at a sim- 
ilar amplitude to the oscillator signal (A mw = 0.0025), 
the intersection occurs halfway between the two signals. 

This is a useful feature that may be of direct bene- 
fit to experimental investigators. The application of mi- 
crowaves to this sort of solid state system is normally via 
a frequency dependent coupling. If multiple frequencies 
are to be used in an experiment, each will need to be 
characterised. This feature will provide an independent 
measure for the amplitude of the microwaves which ac- 
tually couple to the device at different frequencies. It 
should allow the amplitude of the microwaves reaching 
the device to be kept constant even though the coupling 
changes with frequency. 

In addition, it may be possible to characterise the oscil- 
lator in terms of microwave properties, and furthermore 
estimate the effects or significance of the oscillator en- 
ergy transition. Indeed, for high microwave amplitudes 
{A m w > 0.0050), the microwave transition width even- 
tually increases to enclose the oscillator transition, then 
the features merge. 



B. Microwave frequency sweep 

Fig. [7| is presented in a similar manner as Fig. [HJ How- 
ever there are now two frequency axes: the horizontal 
axis represents the frequency of the applied microwave 
drive field, and the vertical axis is the frequency re- 
sponse. It should be remembered that the microwave 
frequency axis is focused near the qubit transition fre- 
quency (f qu bit ~ 3.49GHz) and the diagonally increasing 
line is now the microwave frequency. 

The most interesting features in question are secondary 
splittings, which occur only when the lower Rabi side- 
band passes through the oscillator peak, splitting it. The 
oscillator frequency split is somewhat obscured by the 
large oscillator peak, however the low frequency splitting 
centred on f mw — f osc = 430MHz is clearly visible in the 
magnified section, as indicated by the arrows. 

We find the particular combination required to cause 
maximum splitting only occurs when f Rabl = f mw - f osc 
and the Rabi oscillations have maximum amplitude, (i.e. 
when the microwave is in resonance with the bias set- 
ting). Hence the microwave amplitude and a feature in 
frequency space can be related. The traces created by 
these splittings can be observed at a given applied mi- 
crowave frequency, the traces expand and contract de- 
pendant on the proximity to the splitting. Tracking the 
expansion of these traces it may be possible to tune the 
microwave amplitude for the desired Rabi frequency. 

This could be applied to calibrate the microwave 
waveguide, in which one could determine the effective mi- 
crowave applied to the qubit in terms of the microwave 
amplitude applied at the external end of the waveguide. 
This is important as the characteristics of the waveguide, 
such as impedance, are frequency dependant. If the sep- 
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Oscillator PSD as a function of bias = 0.0025 Oscillator PSD as a function of bias = 0.0050 
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Bias, (2e) Bias, (2e) 



FIG. 6: (Color online) Oscillator PSD as a function of bias, for microwave amplitudes A mw — 0.0025 (A) and A mw = 0.0050 
(B). The red lines track the positions (in frequency) of significant power spectrum peaks (+10dB to +15dB above background), 
the overlaid black and blue lines are the qubit energy and microwave transition (Fig. |3J ■ Unlike Fig. |1J in these figures the 
3.06GHz oscillator circuit can now drive the qubit (Fig. |3J and so creates excitations which mix with the microwave driven 
excitations creating a secondary splitting centred on f mw — f osc (430MHz) . This feature contains the Rabi frequency information 
in the sidebands of the splitting, but now in a different and controllable frequency regime. In addition, the intersection of the 
two differently driven excitations (illustrated in the magnified sections), opens the possibility of calibrating the biased qubit 
against a fixed engineered oscillator circuit, using a single point feature. (s=5x 10 _J ). 



aration of the two known frequencies is fixed, (the oscil- 
lator and the microwave), we know by maximising the 
secondary splitting using the microwave amplitude that 
the Rabi side-band must now be equal to the well defined 
oscillator-microwave frequency separation, and therefore 
the external microwave amplitude that caused it. In ad- 
dition, the maximum splitting solution should indicate a 
correctly biased system. If the oscillator were of a much 
lower frequency (1.36GHz) then only the Rabi oscilla- 
tions and Rabi side-bands would appear (Fig. @}, and 
none of these additional interactions would exist. The 
extra features caused by this small detuning is an impor- 
tant point to reinforce as it links the microwave, oscilla- 
tor, Rabi oscillations and bias field together. 



VI. CONCLUSION 

In this paper we have presented features that exist 
within the power spectra of a charge qubit. The qubit 
is coupled via a single Josephson junction to a (fully 
quantised) oscillator circuit, which provides a bias volt- 
age and a means of measurement. A quantum trajec- 
tory approach has been used, with the quantum state 
diffusion unravelling of the master equation simulating a 
dissipative oscillator, and a quantum jumps unravelling 
simulates spontaneous emissions from the qubit to a lossy 
reservoir. 

A backreaction technique is tested with a quantised 
oscillator to allow the energy level separation of a charge 
qubit to be measured in-situ and with the minimum of 
measuring circuitry and without a separate qubit readout 
device. However care should be taken to avoid detecting 
the higher oscillator harmonics. We also report addi- 
tional features that are created in the power spectrum of 
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FIG. 7: (Color online) Oscillator PSD as a function of the applied microwave drive frequency fmw, for microwave amplitudes 
Amw = 0.0050 (A) and A mw — 0.0100 (B). It is important to notice that there are now two frequency axes per plot, a drive (H) 
and a response (V). Of particular interest is the magnified section which shows clearly the distinct secondary splitting in the 
sub-GHz regime. This occurs due to a high frequency interaction seen in the upper plots, where the lower Rabi sideband of the 
microwave drive passes through the high frequency oscillator signal. The maximum splitting occurs when the Rabi amplitude 
is a maximum, hence this is observed for a very particular combination of bias and drive, which is beneficial for charactering 
the qubit. Most importantly, this would not be observed with a conventional low frequency oscillator configuration as the 
fmw — fosc separation would be too large for the Rabi frequency, (k = 5 x 10~ 5 ). 



the qubit (and the coupled oscillator) when the oscillator 
resonant frequency and microwave frequency are similar 
(3.06GHz and 3.49GHz respectively). These features are 
predominantly discontinuous frequency splittings akin to 
the Autler-Townes effect in quantum optics, and the in- 
tersections of frequency peaks as functions of the control 
fields. These effects are also present at low (sub GHz) 
frequencies which are expected to be more easily accessi- 
ble by experiment than the original 2 to 4 GHz frequency 
range. This frequency range can be moved into an exper- 
imental bandwidth by adjusting the oscillator frequency. 

The close proximity of a second drive field (the 
oscillator), to the microwave drive adds information to 
the noise spectra and allows access to other features, by 
searching for frequency separations and discontinuities 
we can identify bias settings and new amplitudes from 



artefacts in the frequency domain. 
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